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HEAVY QUARK PHYSICS FROM LATTICE QCD 



J.M. Flynn and C.T. Sachrajda 
Department of Physics and Astronomy, University of Southampton 
Southampton SOU IB J, UK 

We review the application of lattice QCD to the phenomenology of b- and e-quarks. 
After a short discussion of the lattice techniques used to evaluate hadronic matrix 
elements and the corresponding systematic uncertainties, we summarise results 
for leptonic decay constants, B—B mixing, semileptonic and rare radiative decays. 
A discussion of the determination of heavy quark effective theory parameters is 
followed by an explanation of the difficulty in applying lattice methods to exclusive 
nonleptonic decays. 



1 Introduction 

The lattice formulation of Quantum Chromodynamics (QCD) and large scale 
numerical simulations on parallel supercomputers are enabling theorists to 
evaluate the long-distance strong interaction effects in weak decays of hadrons 
in general, and of those containing a heavy (c- or b-) quark in particular. The 
difficulty in quantifying these non-perturbative QCD effects is the principal sys- 
tematic error in the determination of the elements of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix from experimental measurements, and in using the 
data for tests of the standard model of particle physics and in searches for 
"new physics" . In this article we will review the main results obtained from 
lattice simulations in recent years, and discuss the principal sources of uncer- 
tainty and prospects for future improvements. More details can be found, for 
example, in the review talks presented at the annual symposia on lattice field 
theory (see Refs.LUlIil) and in references therein. 

In the lattice formulation of quantum field theory space-time is approxi- 
mated by a discrete "lattice" of points and the physical quantities of interest 
are evaluated numerically by computing the corresponding functional inte- 
grals. For these computations to make sense it is necessary for the lattice to 
be sufficiently large to accommodate the particle(s) being studied (L ^> 1 fin 
say, where L is the spatial length of the lattice), and for the spacing between 
neighbouring points (a) to be sufficiently small so that perturbation theory can 
be used to interpolate between the lattice points (oAqcd C 1). The number of 
lattice points in a simulation is limited by the available computing resources; 
current simulations are performed with about 16-20 points in each spatial di- 
rection (up to about 64 points if the effects of quark loops are neglected, i.e. 
in the so called "quenched" approximation). Thus it is possible to work on 
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lattices which have a spatial extent of about 2 fm and a lattice spacing of 0.1 
fm, perhaps satisying the above requirements. 

Since the Compton wavelength of heavy quarks is small, one has to be 
careful in their simulation. With a discretization of the Dirac action, Q(x)(lp + 
mQ)Q(x), where Q represents the field of the quark and ttiq its mass, the 
condition oAqcd <C 1 is not sufficient, one also requires arriQ <C 1. Currently 
available lattices, in simulations performed in the quenched approximation, 
have spacings in the range 1.5 GeV < a -1 < 4 GeV, and so it is not possible to 
simulate 6-quarks directly in this way. Hence, one approach to B-physics is to 
simulate hadrons with heavy quarks which are still considered to be acceptably 
light (typically with masses in the range of that of the charmed quark), and 
then to extrapolate the results to the 6-region, using scaling laws from the 
Heavy Quark Effective Theory (HQET) where applicable. 

An important alternative approach is to use the HQET, and to compute 
physical quantities in terms of an expansion in inverse powers of the mass of 
the 6-quark (and c-quark where appropriate). The first term in this expansion 
corresponds to a calculation using static (infinitely) heavy quarks. The evalu- 
ation of the coefficients in this expansion does not require propagating heavy 
quarks, and hence there are no errors of O(rriQa). However, as will be dis- 
cussed in more detail in Section g the evaluation of the higher order terms in 
this expansion requires the subtraction of power divergences (i.e. terms which 
diverge as inverse powers of the lattice spacing a), which significantly reduces 
the precision which can be obtained. 

In this article we focus on the decays of b- and c-quarks. Another important 
area of investigation is the spectroscopy of heavy quarkoniaJibr a review of 
recent results and references to the originaliiterature see Ref.lia), in which the 
non-relativistic formulation of lattice QCDLLia is generally used. 

The plan of this article is as follows. In the remainder of this Section we 
briefly review the general procedure used for evaluating hadronic matrix ele- 
ments in lattice simulations and discuss the principal sources of uncertainty. 
The following Sections contain the status of the results for the leptonic decay 
constants fg and fjj (Section ||), the .Be-parameter of B°-B° mixing (Sec- 
tion ||), semileptonic decay amplitudes of B- and D- mesons (Section ^) and 
the parameters A (the binding energy of the heavy quark in a hadron), Ai (its 
kinetic energy) and A2 (the matrix element of the chromomagnetic operator) of 
the HQET (Section ||). Section || also contains a general discussion of the the- 
oretical difficulties present in calculations of power corrections, i.e. corrections 
of 0(1/toq), 0(1/toq) etc., in physical quantities. Section || contains a brief 
summary of the status of lattice calculations of exclusive non-leptonic decay 
amplitudes, which is an area of investigation requiring considerable theoretical 
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Table 1: Where in this review to find results for some important physical quantities. 



Quantity 


See. . . 




Leptonic decay constants: fn, /d s , /b, /b s 


Eqs. (pf- 




Renormalization-group-invariant 

hi—TT\£\ T*£) TY1 ptpr TY"1T* r\ — ri TT1 1 Y1Y1 O' ri n 


Eq. © 




z = fB.&£ifB t 8 l £ 


Eq. © 




Form factors /+(0), V(0), A x (0) and A 2 (0) 
for semileptonic D — > if, -K"*,7T, p decays 


Table | 




Form factor / + (0) and decay rate for 
B° -► 7T+r^ 


Table 


(UKQCD values) 


Form factors y(O), Ai(0) and A 2 (0) and rate 
for 5° -> p+r^ 


Table ^ 


(UKQCD values) 


Form factor T(0) = Ti(0) = iT 2 (0) for 


Table || 
and Eq. ( 


(UKQCD values) 



development. Finally in Section |j we present a brief summary and outlook. 
Table [I] shows where our preferred estimates for some of the more important 
physical quantities can be found. 

1.1 Evaluation of Hadronic Matrix Elements 

By using the Operator Product Expansion, it is generally possible to separate 
the short- and long-distance contributions to weak decay amplitudes into Wil- 
son coefficient functions and operator matrix elements respectively. Thus in 
order to evaluate the non-perturbative QCD effects it is necessary to compute 
the matrix elements of local composite operators. This is achieved in lattice 
simulations, by the direct computation of correlation functions of multi-local 
operators composed of quark and gluon fields (in Euclidean space) : 

(0|O(ari,xa,...,aB„)|0)= ^ ([DA^D^D^] e~ s 0(x u x 2 , ... ,x n ) , (1) 
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where Z is the partition function 



Z = J [DA^][D^][Di,]e- s , (2) 

S is the action and the integrals are over quark and gluon fields at each space- 
time point. In Eq. (JjJ) 0(xi, x%, . . . , x n ) is a multi-local operator; the choice of 
O governs the physics which can be studied. 

We now consider the two most frequently encountered cases, for which 
n=2 or 3. Let 0{x\,X2) be the bilocal operator 

02(11,3:2) - T{J h ( Xl )jl(x 2 )} , (3) 

where Jh is an interpolating operator for the hadron h whose properties we 
wish to study. In lattice computations we evaluate the two point correlation 
function 

C 2 (U = £<0|O2(*.0)|0> ■ (4) 

X 

For sufficiently large positive t x one obtains: 

p—mutx 

C 2 (t x )~- |(0|J fc (0)|ft)| 2 . (5) 

2uih 

In Eq. (Q) mh is the mass of the hadron h, which is assumed to be the lightest 
one which can be created by the operator J^. The contribution from each 
heavier hadron, h! with mass my say, is suppressed by an exponential factor, 
exp ( — (rrih' — rrih)tx)- In lattice simulations the correlation function Ci is 
computed numerically, and by fitting the results to the expression in Eq. (||) 
both the mass to/j and the matrix element (0 | J/i(0)| h) can be determined. In 
this case the hadron h is at rest, but of course it is also possible to give it a 
non-zero momentum, p say, by taking the Fourier transform in Eq. (^) with 
the appropriate weighting factor exp(ip-x). 

As an example consider the case in which h is the i?-meson and Jh is the 
axial current (with the flavour quantum numbers of the B- meson) . In this 
case one obtains the value of the leptonic decay constant /g, 

(0|^(0)| J B(p)) = / B p /t . (6) 

It will also be useful to consider three-point correlation functions: 

C 3 (t x ,t v ) =5> 4 P- X e^(0| J 2 (x,t x )O(y,^)J 1 t (0,0)|0) , (7) 
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where, J\ and J 2 are the interpolating operators for hadrons h\ and /i 2 respec- 
tively, O is a local operator, and we have assumed that t x > t y > 0. Inserting 
complete sets of states between the operators in Eq. (0) we obtain 

Cat*,,*,) = -^r ^ (0|J 2 (0 ) 0)|ft 2 (p, J E 2 ))x 

</i 2 (p, £ 2 )|O(0, 0)|/n(p+q, Ex)) (fti(p+q, jt(0, 0)|0> + • • ■ , (8) 



where = y/ m 2 + (p+q) 2 , -B 2 = \/ to 2 , + p 2 and the ellipsis represents the 
contributions from heavier states. The exponential factors, exp(— E\t v ) and 
exp ( — E 2 (t x — t y )), assure that for large time separations, t y and t x — t y , the 
contributions from the lightest states dominate. All the elements on the right- 
hand side of Eq. (||) can be determined from two-point correlation functions, 
with the exception of the matrix element (/i 2 |0|/ii). Thus by computing two- 
and three-point correlation functions the matrix element (/i 2 |0|/ii) can be 
determined. 

The computation of three-point correlation functions is useful in studying 
semileptonic and radiative weak decays of hadrons, e.g. if h\ is a B-meson, hi 
a D meson and O the vector current Jyy^c, then from this correlation function 
we obtain the form factors relevant for semileptonic B — ► D decays. 

We end this brief summary of lattice computations of hadronic matrix 
elements with a word about the determination of the lattice spacing a. It 
is conventional to introduce the parameter j3 — 6/gQ (a), where go(a) is the 
bare coupling constant in the theory with the lattice regularization. It is (3 
(or equivalently (70(a)) which is the input parameter in the simulation, and 
the corresponding lattice spacing is then determined by requiring that some 
physical quantity (which is computed in lattice units) is equal to the exper- 
imental value. For example, one may compute m p a, where m p is the mass 
of the p-meson, and determine the lattice spacing a by dividing the result by 
769 MeV. 



1.2 Sources of Uncertainty in Lattice Computations 

Although lattice computations provide the opportunity, in principle, to eval- 
uate the non-perturbative QCD effects in weak decays of heavy quarks from 
first principles and with no model assumptions or free parameters, in practice 
the precision of the results is limited by the available computing resources. In 
this section we outline the main sources of uncertainty in these computations: 

• Statistical Errors: The functional integrals in Eq. (Q) are evaluated by 
Monte-Carlo techniques. This leads to sampling errors, which decrease as 
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the number of field configurations included in the estimate of the integrals 
is increased. 

• Discretization Errors: These are artefacts due to the finiteness of the lat- 
tice spacing. Much effort is being devoted to reducing these errors either 
by performing simulations at several values of the lattice spacing and ex- 
trapolating the results to the continuum limit (a = 0), or by "improving" 
the lattice formulation of QCD so .thai, the error in a simulation at a 
fixed value of a is formally smaller EinI3. In particular, it has recently 
been shown to be possible to formulate lattice QCD in such a way that 
the discretization errors vanish-quadratically with the lattice spacing LLa, 
even for non-zero quark massesE3. This is in contrast with the traditional 
Wilson formulation EJ, in which the errors vanish only linearly. Some of 
the simulations discussed in the following sections have used a tree-level 
improved action and operators; in these studies the artefacts formally 
vanish more quickly (like aa s (a)) than for the Wilson action. In the fol- 
lowing, this tree-level improved action will bfi_denoted as the SW (after 
Shcikholcslami-Wohlert who first proposed itll3) or "clover" action. 

• Extrapolations to Physical Quark Masses: It is generally not possible to 
use physical quark masses in simulations. For the light (u- and d-) quarks 
the masses must be chosen such that the corresponding 7r-mesons are suf- 
ficiently heavy to be insensitive to the finite volume of the lattice. In 
addition, as the masses of the quarks are reduced the computer time nec- 
essary to maintain the required level of precision increases rapidly. For 
the heavy quarks Q (i.e. for c, and particularly for b) the masses must be 
sufficiently small so that the discretization errors, which are of 0(mQa) 
or 0(niQa 2 ), are small. The results obtained in the simulations must then 
be extrapolated to those corresponding to physical quark masses. 

• Finite Volume Effects: We require that the results we obtain be inde- 
pendent of the size of the lattice. Thus, in principle, the spatial size of 
the lattice L should be ^ I fm (in practice L > 2-3 fm) , and as men- 
tioned above, we cannot use very light u- and d-quarks (in order to avoid 
very light pions whose correlation lengths, i.e. the inverses of their masses, 
would be of O(L) or greater). 

• Contamination from Excited States: These are the uncertainties due to 
the effects of the excited states, represented by the ellipsis in Eq. (||). In 
most simulations, this is not a serious source of error (it is, however, more 
significant in computations with static quarks). Indeed, by evaluating 
correlation functions (Jh{x) ^(0)) for a variety of interpolating operators 
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Figure 1: Comparison of lattice spacings (the quantity plotted is actually the inverse lattice 
spacing, a -1 ) determined from different physical pmantities in quenched and unquenched 
simulations by the SESAM and TxL collaborationscj- F denotes the value determined from 
the static quark potential and p denotes the value determined from the p-meson mass, while 
2S-1S and IP— IS 1 denote values obtained from energy level splittings in QQ bound states 
using the lattice formulation of nonrelativistic QCD. 



{Jh, J'h}, it is possible to obtain the masses and matrix elements of the 
excited hadrons (for a recent example seecj). 

• Lattice- Continuum Matching: The operators used in lattice simulations 
are bare operators defined with the lattice spacing as the ultra-violet cut- 
off. From the matrix elements of the bare lattice composite operators, 
we need to obtain those of the corresponding renormalized operators de- 
fined in some standard continuum renormalization scheme, such as the 
MS scheme. The relation between the matrix elements of lattice and 
continuum composite operators involves only short-distance physics, and 
hence can be obtained in perturbation theory. Experience has taught us, 
however, that the coefficients in lattice perturbation theory can be large, 
leading to significant uncertainties (frequently of 0(10%) or more). For 
this reason, non-perturbative techniques to evaluate the renormalization 
constants which relate the bare lattice operators to the renormalized ones 
have been developed, using chiral Ward identities where possibJeEj or by 
imposing an explicit renormalization conditions (see also Refs.E3o), thus 
effectively removing this source of uncertainty in many important cases. 

• "Quenching" : In most current simulations the matrix elements are evalu- 
ated in the "quenched" approximation, in which the effects of virtual quark 
loops are neglected. For each gluon configuration {[/^(x)}, the functional 
integral over the quark fields in Eq. (Q) can be performed formally, giv- 
ing the determinant of the Dirac operator in the gluon background field 
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Figure 2: Diagram representing the leptonic decay of the B-meson. 

corresponding to this configuration. The numerical evaluation of this de- 
terminant is possible, but is computationally very expensive, and for this 
reason the determinant is frequently set equal to its average value, which 
is equivalent to neglecting virtual quark loops. Gradually, however, un- 
quenched calculations are beginning to be performed, e.g. in Figure |l] we 
show the lattice spacing obtained by the SESAM and T%L collaborations 
from rfour physical quantities in both quenched and unquenched simula- 
tionscJ. In the quenched case there is a spread of results of about ±10%, 
whereas in the unquenched case the spread is smaller (although the errors 
are still sizeable for some of the quantities used). In the next 3-5 years it 
should be possible to compute most of the physical quantities discussed 
below without imposing this approximation. 

2 The Leptonic Decay Constants fr> and fs 

In this section we review the current status of calculations of the leptonic de- 
cay constants fjj and Jb- Leptonic decays of heavy mesons, see Figure ^|, are 
particularly simple to treat theoretically |]. In each case the strong interac- 
tion effects are contained in a single parameter, called the decay constant, fo 
or fs- Parity symmetry implies that only the axial component of the V—A 
weak current contributes to the decay, and Lorentz invariance that the matrix 
element of the axial current is proportional to the momentum of the meson 
(with the constant of proportionality defined to be the decay constant) as, for 
example, in Eq. (|^). 

Knowledge of fs would allow us to predict the rates for the corresponding 
decays: 



where the 0(a) corrections are also known. The value of fs is very important 

"For simplicity, the presentation here is for pseudoscalar mesons D and B. A parallel 
discussion holds also for the vector mesons D* and B* . 



t{b -> i n + iv a ) 




B 



(l + 0(a)), (9) 
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Table 2: Results for fu and fjj using the conventional formulation of heavy quarks. All 
values were obtained in the quenched approximation, except for those from the HEMCGC 
collaboration which were computed for two flavours of sea quarks (n/ = 2) (we have combined 
quoted systematic errors from HEMCGC in quadrature). Infinite values of /3 symbolise 
results obtained after extrapolation to the continuum limit. The error quoted by APE 
includes an estimate of discretization uncertainties taken from a comparision of results at 
j3 = 6.0 and 6.2. Results from MILC and JLQCD are preliminary. 
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norm 


fa/MeV 
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/d 3 1 Jd 


MILC 


26 


97 
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nr 


186(10)(?g)(8) 


199(8)(f?)( 1 g) 


1.09(2)(f)(g) 


JLQCD 


27 


97 


00 
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nr 
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APE 
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00 





nr 
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nr 
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1.09(2)(5) 


UKQCD 


31 


94 


6.2 
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185(f)( 4 ?) 
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UKQCD 


31 


94 


6.0 


1 


rcl 




225(ig)(l§) 


1.13(f) 


BLS 


:i2 


94 


6.3 





nr 


208(9)(35)(12) 


230(7)(30)(18) 


1.11(6) 


HEMCGC 


33 


94 


5.3 





nr 


215(5)(53) 


287(5)(60) 




HEMCGC 


34 


93 


5.6 







200-287 
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ELC 


35 


92 


6.4 
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210(40) 


230(50) 




PCW 


■M> 


91 


6.0 
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209(18) 




ELC 
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37 
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6.0 
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197(14) 
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DeGL 
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88 


6.0 





rcl 
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222(16) 


1.17(22) 
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:i!t 


88 


6.1 





rcl 


174(26)(46) 


234(46)(55) 


1.35(22) 



in describing B-B mixing as explained in Section |3_|. Knowing fs would also 
be useful for our understanding of other processes in B-physics, particularly 
those for which "factorization" turns out to be a useful approximation. 

Lattice results for decay constants of charmed and bottom mesons obtained 
over the last ten years are summarised in Tables ^| and |j| These are presented 
using a normalization in which f 7T + ~ 131 MeV. 

Since we are dealing with a heavy quark with mass mg, the product 
m,Qa is large and can be an important source of mass-dependent discretization 
errors. This has been addressed in two ways in the decay constant calcu- 
lations from 1994 onwards. Some studies, denoted by c sw = 1, have used 
the Sheikholeslami-Wphlert (SW or clover) improved action E£l together with 
improved operators til, which removes tree level 0(a) errors, leaving ones of 
0(a s iriQa). In the near future we can look forward to the application of re- 
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Table 3: Results for /b and fs s using the conventional formulation of heavy quarks. All 
values were obtained in the quenched approximation, except for those from the HEMCGC 
collaboration which were computed with two flavours of sea quarks (nf = 2) (we have 
combined quoted systematic errors from HEMCGC in quadrature). Infinite values of f) 
symbolise results obtained after extrapolation to the continuum limit. The error quoted by 
APE includes an estimate of discretization uncertainties taken from a comparision of results 
at /3 = 6.0 and 6.2. Results from MILC and JLQCD are preliminary. 
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205(40) 
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cently developed techniques to reduce the discretization errors still further, 
to ones of Oimqa 2 ) The second approach introduces a revised normal- 

ization of the quark fields in simulations using the standard Wilson fermion 
action. This is designed to remove higher order effects in rriQa from the heavy 
quark propagator, but only at tree level in the strong interactions, and is 
distinguished in the table by the label 'nr', denoting a nonrelativistic normal- 
ization (in contrast to the standard 'rel' or relativisitic one). The theoretical 
significance of this normalization factor is not completely clear; in particular, 
it does not eliminate all the 0{m,Qa) effects. The nonrelativistic normaUza±km 
is often denoted 'KLM' in the literature, after some of its originators&c2lc3. 
For more details of the lattice formulations and improvement procedures, see 
the recent review by Wittigci 

The main efforts of the lattice community are being devoted to controlling 
the systematic uncertainties, which dominate the errors. This is illustrated 
by the latest- and still preliminary, results, which come from the MILC o 
and JLQCDo collaborations. These collaborations are carrying out extensive 
simulations using many different lattice spacings to allow them to extrapolate 
to the continuum limit. JLQCD study different prescriptions for reducing the 
0(m,Qa) discretization errors associated with heavy quarks with the aim of 
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demonstrating that all results converge in the continuum limit. Their error is 
statistical combined with the uncertainty due to the spread over prescriptions. 
In Figure ^ we show the continuum extrapolation of their results for fu and 
/b for two different lattice fermion formulations. In this case the continuum 
values of fs and fo obtained using the two formulations agree remarkably well; 
the agreement is still acceptable (although not so remarkable) when quantities 
other than the string tension are used to determine the lattice spacing. We also 
note that in this case the dependence on the lattice spacing is milder for the 
improved action as expected (although further studies are needed to confirm 
whether this is a general feature). 

MILC simulate a range of heavy quark masses, giving meson masses strad- 
dling rriD, together with a static (infinite mass) quark, allowing an interpola- 
tion to tub- The continuum extrapolation of their results is illustrated in 
Figure In the MILC results, the first error is statistical and the second 
is systematic within the quenched approximation. MILC also have some un- 
quenched simulation results (with nj = 2 dynamical flavours). Their final 
error is for the effects of quenching, which they estimate from: (i) the differ- 
ence between the unquenched results at the smallest available lattice spacing 
and the quenched results interpolated to the same spacing, and (ii) comparing 
the results when the lattice spacing is fixed by /„■ or m p .| Their results sug- 
gest that unquenching may raise the value of the decay constants. This agrees 
with estimates c£la, using the difference between chiral loop contributions in 
quenched and unquenched QCD, that fs and fs 3 in full QCD are increased 
by 10-20% over their quenched values. Lattice simulations are possible using 
scalar fields with a quark-like action which act effectively as negative numbers 
of ordinary quark flavoursLJ. Calculations of Jb in the static limit extrapolated 
from negative numbers of flavours also suggest an increase of about 20%. 

At this point we feel we need to add a word of caution. In their evalua- 
tion of /s, both the JLQCD and MILC collaborations simulate heavy quarks 
with masses uiq satisfying niQa > 1 and apply a nonrelativistic or 'KLM' nor- 
malization prescription. Even with such a normalization, however, one should 
worry about the discretization errors when such large masses are used. Not 
only would one expect these errors to be sizeable at a fixed value of beta, but 
the extrapolation to the continuum limit is now complicated since all powers of 
rriQa contribute significantly to the results which are being extrapolated. Set 
against these worries is the observation by the MILC collaboration that they 
obtain very similar values for fs if they restrict the set of points to be extrapo- 
lated to those with heavy quarks with reasonably small masses (the differences 

'For / n and fg they also compare results where the strange quark mass is determined 
by m K or m^,. 
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Figure 3: Continuum extrapolation of JLQCD results (preliminary) for fr> and fs- The 
graph shows agreement between results from two different lattice formulations when extrap- 
olated to zero lattice spacing. Open symbols denote an unimproved Wilson action (c sw = 0), 
filled symbols denote an improved Sheikholeslami-Wohlert or Clover (c sw = 1) action. 
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Figure 4: Continuum extrapolation of MILCcj results (preliminary) for fg. Circles (tri- 
angles) denote results from quenched (unquenched nj = 2) simulations. The solid line and 
square denote the continuum extrapolation by a linear fit in the lattice spacing a to the 
quenched results, while the dashed line and diamond denote a weighted average of the three 
quenched results at smallest a. The scale is set by f w . 
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are included in their error bars). Our own strong preference is to restrict the 
simulations to heavy quarks with masses significantly smaller than 1 in lattice 
units (and to use improved actions). 

In obtaining a continuum value for Jb one has to decide in what order to 
perform the extrapolations in a and tuq. The MILC collaboration obtain /g 
by fitting, at each finite lattice spacing a, the value of the decay constant of a 
static quark (infinite mass) and those of several mesons with finite masses (it 
should be noted that these have different systematic errors). They then extrap- 
olate the results for / r to th e continuum limit to obtain their final result. An 
alternative procedure E3oo is to perform separately the continuum extrapo- 
lations of the static results and results obtained at fixed finite meson masses. 
The mass dependence of the continuum points is then fitted to determine the 
physical decay constants. The most recent analysis using this method has been 
performed by WittigE-3, combining lattice data from many simulations. The 
MILC and JLQCD data is excluded, since it is still preliminary, and the results 
obtained with the clover action have also not been included since it was felt 
that there is insufficient data to perform the continuum extrapolation. This 
analysis shows some evidence for discretization errors at large quark masses, 
where results using different normalization procedures, mass definitions and so 
on, begin to diverge. Wittig's final fit, therefore, is to the static point and to 
points, with meson masses less than about 2 GeV. The results of this procedure 
giveEJ 

f D = 191±iMeV (10) 
f Ds = 206+iMeV (11) 
f B = 172±§lMeV (12) 

The statistical error and the systematic error arising from fixing the lattice 
spacing are the dominant (and comparable) uncertainties in these numbers. 
The central values are obtained using the continuum extrapolated value of 
to fix the lattice spacing a. 

The dimensionless ratios fo B / Jd and fg t / fs are also of great interest and 
have the advantage that many systematic effects should partially cancel. The 
lattice results illustrate, as expected, that the value of the decay constant de- 
creases as the mass of the light valence quaik decreases. Here a straightforward 
extrapolation of world results by WittigE.il gives 

f D.I in = 1.08 ±0.08, (13) 
SbJSb = 1.14 ±0.08, (14) 

which are compatible with the preliminary MILC results H These results are 
obtained in the quenched approximation. MILC estimate the quenching error 
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in the ratios to be about 3%. One would expect a small correction here from 
cancellation of the effects in Jb s and fg separately. We note that the chiral 
loop estimate for the difference of the ratio fs s / Jb from 1 can be surprisingly 
largec3o, e.g. for one choice of parameters it can be as large as +16%. 

2.1 Summary and Conclusions 

We now summarise this rather lengthy discussion and present our conclusions. 
The two principle sources of uncertainty in the calculations of the decay con- 
stants are those due to discretization errors and to quenching. The effects of 
the latter are only now beginning to be studied, and early indications suggest 
that the values of the decay constants may increase when the effects of quark 
loops are fully included, but this still needs to be confirmed. The discretization 
errors are currently studied either by using improved actions and operators, or 
by performing the computations at a sequence of lattice spacings and extrapo- 
lating the results to the continuum limit. The results from the two approaches 
are consistent, and further confirmation that these errors are under control 
will come afte r the implementation of the 0(a 2 ) improvement techniques pro- 
posed in Refs.liao. Our view of the current status of the calculations can be 
summarised by the following values for the decay constants and their ratios: 



f D = 200±30MeV (15) 

f Ds = 220±30MeV (16) 

f B = 170±35MeV (17) 

f Bs = 195±35MeV (18) 

IdJId = 1.10 ±0.06 (19) 

SbJSb = 1.14 ± 0.08 (20) 



The principal difficulty in presenting results for physical quantities such as 
those in Eqs. (|l5|-|2C|) is to estimate the errors. The value of the lattice spacing 
typically varies by 10% or so depending on which physical quantity is used to 
calibrate the lattice simulation. This variation is largely due to the use of the 
quenched approximation. We therefore consider that ±10% is an irreducible 
minimum error in decay constants computed in quenched simulations (15% 
when they are computed in the static limit, since in that case it is fsy/rnB 
which is computed directly). The remainder of each error in Eqs. dl5|-p0|) and 
below is based on our estimates of the other uncertainties, particularly those 
due to discretization errors and the normalization of the lattice composite 
operators. As explained in the Introduction, much successful work is currently 
being done to reduce these uncertainties. 
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It is very interesting to compare the lattice prediction for fp 3 in Eq. (16) 
with experimental measurements. Combining four measurements of /c^. from 
D s — > (ii> decays, the rapporteur at the 1996 ICHEP conference foundcj 

f Da = 241 ± 21 ± 30 MeV. (21) 

In spite of the sizeable errors, the agreement with the lattice prediction is very 
pleasing and gives us further confidence in the predictions for fp and related 
quantities. 

We end this section with a comment on the validity of the asymptotic 
scaling law for the decay constants. For sufficiently large masses of the heavy 
quark, the decay constant of a heavy-light pseudoscalar meson P scales with 
its mass mp as follows: 



a s (m P 



-^c {l + 0Mmp))}+0 Q_) 



(22) 



where A is independent of mp. Using the leading term of this scaling law, a 
value of 200 MeV for fp would correspond to fp — 120 MeV. Results from 
the lattice computations shown in Table ^, however, indicate that fp is signif- 
icantly larger than this and that the 0(1/ mp) corrections on the right-hand 
side of Eq. ( p2] ) are considerable. The coefficient of the 0(l/mp) corrections 
is typically found to be between 0.5 and 1 GeV, on the large side of theoretical 
expectations. 



3 B° - B° Mixing 

B-B mixing provides important constraints for the determination of the angles 
of the unitariiy triangle (for a review and references to the original literature 
see e.g. Ref.Ei). In this process, the strong interaction effects are contained in 
the matrix element of the A£?=2 operator: 

M{ii) = (B°\ 67^(1-75)9 67 M (l-75)g|5°), (23) 

where q—d or s (unless indicated otherwise, we assume that q=d). The argu- 
ment n implies that the operator has been renormalized at the scale \x. It is 
conventional to introduce the £?e-parameter through the definition 

M(») = ^flm B B B (n). (24) 

The dimensionless quantity Bp is better-determined than fp in lattice calcu- 
lations, so that the theoretical uncertainties in the value of the matrix element 
M, needed for phenomenology, are dominated by our ignorance of fp. 
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Table 4: Results for the mixing parameter Bb obtained in the quenched approximation 
using static fo-quarks (the value for c sw refers to the light quark action). Bg(mi,) with 
mi, = 5 GeV is calculated from the raw lattice matrix elements and then converted to Bg° , 
the renormalization group invariant B parameter. For ease of comparison we have used 
the same choice of parameters in all cases together with NLO matching from lattice to 
continuum. Where this has changed the results from those quoted by the authors, the result 
obtained by us is indicated in oblique type. Ml, M2 and M3 denote three different but 
equivalent ways of organising the matching calculations (see text). 





Yr 


a 






B B {m b ) 


B B 


gmEI 


97 


6.0 


1 


Ml 
M2 
M3 


0.76(5) 
0.54(4) 
0.86(5) 


1.21(8) 
0.86(6) 
1.37(8) 


ukqcdEI 


96 


6.2 


1 


Ml 
M2 
M3 


0.76(6) 
0.57(6) 
0.82(6) 


1.20(9) 
0.91(9) 
1.31(9) 


KenH 


96 


6.0 





Ml 
M2 
M3 


0.93(4) 
0.75(4) 
0.98(4) 


1.48(7) 
1.19(7) 
1.56(7) 



Bs(fJ-) is a scale dependent quantity for which lattice results are most often 
quoted after translation to the MS scheme. The next-to-leading order (NLO) 
renormalization group invariant .B-parameter (Bg°) is defined by 

Bt = *s(»r 2/0 ° (l + ^Uv) B B M, (25) 

where [3o — 11 — 2n//3 and J n . is obtained_£com the one- and two-loop anoma- 
lous dimensions of the AB—2 operator by Ell, 

J "^^( 7 °l~ 71 )' (26) 

with j3\ — 102 — 38rt//3, 70 = 4 and 71 = — 7 + 4n//9. In the discussion below 
we choose [i — m&. 

Results for Bb are summarised in Table || for static b quarks and in Table || 
for propagating 6-quarks. We now discuss these two cases in turn. 

3.1 Bb obtained using static b-quarks 

The main difficulty in trying to determine the value of Bb using static heavy 
quarks is due to the large perturbative corrections which are encountered when 
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calculating the matrix element M(fi) in some standard renormalization scheme 
(e.g. the MS scheme) in full QCD from those measured on the lattice in the 
effective theory. We will illustrate this below. Because of the relatively large 
uncertainties which result from the truncation of the perturbation series, we 
believe that, at present, the calculations of Bb in the static theory add little to 
the information obtained with propagating quarks. For this reason our "best" 
lattice value for Bb is obtained using results obtained with propagating quarks 



only, and will be presented in section 3.2 below 



In order to compute M(mb) one needs to evaluate the matrix elements 
of four AB = 2 operators in the static lattice theory (denoted by Oi(a) with 
i = L, N, R, S). The relation between M(m&) and the matrix elements of the 
Oi is of the form: 

M(m b ) = 

Zi(rrib, a) (Oi(a)) . (27) 

i=L,N,R,S 

The explicit form of the operators and details of the matching can be found in 
Refs.E3T3. The matching can be performed in three stages: 

(i) From the matrix elements of the lattice operators in the effective theory 
one can obtain the corresponding ones renormalizcd in some standard 
continuum scheme at a scale n = a^ 1 . 

(ii) Using the renormalization group equations one then obtains the contin- 
uum effective theory matrix elements of operators renormalized at a scale 
/i = m b . 

(iii) Finally, by matching the continuum effective theory and full QCD one 
obtains M(m b ). 

We illustrate this procedure by considering Zl, which in practice is the most 
complicated case (the operator Ol is just that in Eq. ( p^ ) but with static heavy 
quarks). It can be written in the form: 

Oi s {m b )\ dl f a s (a _1 ) - a s (m b ) 



Cs(m b ) 



a s (mb)\ 2 ( a s (m b ) \ dl I \ ( a s (a 1 



a s (a l ) J \a s (a 



K}[l + ^^D L ) ,(28) 



where J and K are constants obtained from the anomalous dimension matrix 
of the operators Ol and Os in the continuum effective theory, and Dl is 
a constant (which arises in step (i) of the matching procedure). The terms 
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containing both the scales mj and a -1 in Eq. ( p8f ) arise from the evolution in 
step (ii) and the coefficients Cl and Cs from the matching in step (iii). 

To obtain Bb we divide the result for M(rrib) by f B , and so in the denom- 
inator there is a factor of Z A , where Za is the constant relating the physical 
axial current to that in the lattice static theory. Za is also calculated in per- 
turbation theory. 

We consider three methods which have recently been used to treat the 
perturbative corrections; these methods differ only at higher orders in per- 
turbation theory for which the coefficients are unknown. They are therefore 
equivalent at NLO. Because some of the coefficients in the Z^s are large, this 
leads to significant uncertainties in the result for Bb- The three methods are: 

Ml: the expressions for each of the Zi are evaluated exactly as in the example 
ofEq. d|), 

M2: the terms of 0{a 2 s ) in each Zi are dropped (i = L, N, R, S, A), 

M3: the terms of 0{a 2 s ) in each Zi/Z\ are dropped (i — L, N, R, S). 

Ml and M2 are the methods used by Gimenez and Martmeili in El, while M3 
is the "fully-linearised" method of the Kentucky groupEj. 

For the numerical estimates in Table f| we adopt the same choice of pa- 
rameters as Gimenez and Martinelli inE3, using A^L = 200 MeV, nt = 4 and 

' MS „ 1 

mb = 5 GeV for the lattice-to-continuum matching to find Bsimb) [|. To obtain 
B B l ° we set rif = 5 in Eq. ( |25| ) and keep a s continuous at the threshold at nib- 
Where necessary, we have taken the raw matrix elements from the different 
groups and have simply propagated the errors in the linear combination which 
determines Bb- 

As Table |] shows, because of the large coefficients, there are significant 
differences between the formally equivalent procedures Ml, M2 andM3. Fur- 
thermore, the results from UKQCDcil and Gimenez and Martinelli t3, which 
use the same action (i.e. the tree-level improved SW action corresponding to 
c sw — 1) for the light quarks, agree when subjected to the same analysis pro- 
cedure, whereas the Kentucky group resultsE3, using the unimproved Wilson 
action (c sw = 0) but with tadpole-improved perturbation theory in the match- 
ing, are significantly larger. A similar conclusion is reached by Wittig inca. 
Because of these uncertainties, as mentioned above, we will not make further 
use of these results, but will use those obtained with propagating quarks to 
obtain the value of Bb- 

c In step (i) of the maiiehing procedure we use cty(q*a) evaluated at q* = 2.18a — 1 as the 
strong coupling constant e3. 



19 



Table 5: Results for the mixing parameter Bb obtained in the quenched approximation using 
propagating 6-quarks. The authors' results for Bg(fi) at scale fj, or for Bg(mi,) with = 
5GeV are quoted. They are then scaled to tnj and converted to Bj°, the renormalization 
group invariant B parameter. The results may differ slightly from those in the original 
articles owing to the particular choice of parameters used here: authors' numbers are shown 
in upright type, numbers produced by us in oblique type. 







Yr 




fj,/GeV 


Bb{h) 


B B (m b ) 


DIllO 

U B 


BBS 


30 


97 


00 


2 


1.02(13) 


0.96(12) 


1.53(19) 


JLQCD 


31 


96 


6.3 






0.840(60) 


1.34(10) 


JLQCD 


31 


96 


6.1 






0.895(47) 


1.42(7) 


B£ 


1 


96 


00 


2 


0.96(6)(4) 


0.90(6) (4) 


1.44(9)(6) 


ELC 


i5 


92 


6.4 


3.7 


0.86(5) 


0.84(5) 


1.34(8) 


BDHS 


if) 


88 


6.1 


2 


1.01(15) 


0.95(14) 


1.51(22) 



3.2 Bb obtained using propagating b-quarks 

Calculations with propagating heavy quarks are reported in Table H. We have 
scaled all results to a common scale mt = 5 GeV, with the same parameter 
choice as for the static results above, and then converted both to Bb at NLO. 
The results show np_pbservable dependence on the lattice spacing, although 
the authors of Ref. perform a linear extrapolation to the continuum limit, 
which is the reason for the relatively large error in the corresponding result. 
Although the quoted errors are largely statistical, the different groups treat the 
perturbative matching in different ways, so that it is not appropriate to simply 
perform a weighted average. Our preferred estimate based on the results in 
Table | is 

B B ° = 1.4(1) . (29) 

In estimating the error we have assumed that the results are almost inde- 
pendent of the lattice spacing and have not tried to quantify the effects of 
quenching fl. Calculations to be performed in the near future will use improved 
actions and operators (reducing the discretization errors to ones which are 
quadratic in the lattice spacing) and will avoid lattice perturbation theory by 
using non-perturbative renormalization. 

The relevant quantity for B-B mixing is f B Bs- Taking the result in 
Eq. (H) above for B B la with f B = 170(35) MeV from Eq. (0) gives 

f B = 201 (42) MeV (30) 

d The second error in the result from Ref. [] is the authors' estimate of the quenching 
errors. 
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as our lattice estimate. An interesting dimensionless quantity is the ratio 

f Bd VB Bd 



For propagating quarks T jcpmbining the result f B Jf B = 1-14(8) from Eq. (20) 
with BbJBb = 1.00(3)E3 gives 

£ = 1.14(8). (32) 

A recent direct extraction af-jthe matrix element M(/j.) gives the ratio r s d ^ 
£ 2 m| /m Bd = 1.54(13) (32) E3. In the static case, Gimenez and Martinelli E3 
find r sd = 1.43(7) by combining f B Jf B = 1.17(3) and B B JB B = 1.01(1), 
and r s d = 1.35(5) from a direct evaluation of the four-quark matrix element 
measured on the same gauge configurations. The results of the two methods 
are quite consistent, but future calculations should improve on the precision 
of r sd . 

The above results are largely from quenched calculations. For B Bs /B B , 
numerical evidence suggests a small increase on two-flavour dynamical config- 



urations □ but the chiral loop estimate EJe3 is for a decrease of —0.04 in the 
ratio. We must wait for reliable simulations with dynamical quarks before the 
size of quenching effects can be determined with confidence. 

4 Semileptonic Decays of D and .B-mesons 

We now discuss semileptonic decays of D and B-mesons, considering in turn 
the cases in which the c-quark decays into an s- or ci-quark and the 6-quark 
decays into a c- or it-quark. The B-meson decay is represented in Figure gj, 
though the diagram could simply be relabelled to describe the D-meson decay. 
It is convenient to use space-time symmetries to express the matrix elements 
in terms of invariant form factors (using the helicity basis for these as defined 
below). When the final state is a pseudoscalar meson P, parity implies that 
only the vector component of the V—A weak current contributes to the decay, 
and there are two independent form factors, / + and / , defined by 



(P(k)\V»\B(p)) = f + (q 2 



+ /V) m| 2 m V , (33) 

where q is the momentum transfer, q — p—k, and B{p) denotes either a B or 
D meson. When the final-state hadron is a vector meson V, there are four 
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V-A 




D, D*, 7v, p 



Figure 5: Diagram representing the semileptonic decay of the B-meson. q represents the 
light valence antiquark, and the black circle represents the insertion of the V—A current with 
the appropriate flavour quantum numbers. 



independent form factors: 

(V(k,e)\V»\B(p)) = Wl f ] e^efakg (34) 
tub + my 

(V(k,e)\A^\B(p)) = i(m B +m v )A 1 (q 2 )e*^ -i A ^ e*-p{p+ky 

mB+mv 

+ i^J-2m V E*-pq> i , (35) 
T 

where e is the polarization vector of the final-state meson, and q = p—k. Below 
we shall also discuss the form factor A , which is given in terms of those defined 
above by Aq = A + A3 , with 

mB+mv . m B -m v 

A3 = — ^ Ai A 2 . 36 

Zmv Lrtiv 



4-1 Semileptonic D Decays 

The decays D — ► Kl + vi and D — > K*l + vi provide a good test for lattice calcu- 
lations since the relevant CKM matrix element V cs is well constrained in the 
standard model. The form factors for the decays D — > nl + vi and D — > pl + vi 
can also be computed in lattice simulations. As explained in the introduction, 
charm quarks are light enough to be simulated directly (though one still needs 
to be wary of mass-dependent discretization errors). Furthermore, strange 
quarks can also be simulated directly, so for D — » K or K* decays there is 
only one quark for which a chiral extrapolation needs to be performed. For 
semileptonic D-meson decays the whole physical phase space can be sampled^, 

e ln addition, one obtains the form factors for unphysical, negative, values of q 2 . 
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-0.1 0.0 0.1 0.2 

Figp*e 6: UKQCD lattice results for the form factor f + (q 2 ) in semileptonic D — > K de- 
cayE3- The form factor is plotted as a function of the dimensionless variable q 2 a 2 and before 
multiplication by the appropriate lattice-to-continuum normalization factor Z^F . The label 
K l = K crit! K ! s = K s shows that the light quark masses have been extrapolated to give the 
physical D — + K form factor. The curves show fits of pole dominance behaviour to the data: 
the solid line is a two parameter fit, the dashed line a fit where the pole mass is fixed to the 
corresponding 1 _ vector meson mass determined in the same lattice calculation. The cross 
and diamond show the interpolations of the fits to q 2 = 0. 

while keeping the spatial momenta of the initial and final state mesons small 
in order to minimise the momentum-dependent discretization errors. 

Although the lattice calculations actually measure the q 2 dependence of 
the fornipiactors, as shown by the example of / + for the D — > K decay in 
Figure we follow the standard practice of quoting values at q 2 = 0. In 
contrast to the case for B decays to be discussed below, we emphasise that 
this involves an interpolation and so is relatively well controlled. We refer the 
reader to the original papers for detailed discussions of the q 2 -dependence of 
the form factors. 

Lattice results for the D — > form factors are collected in Table |^ and 
illustrated in Figure [?], while results for the D — > 7T, p form factors appear in 
Table ^. These are all from quenched simulations and no group has performed 
a continuum extrapolation. Our summary view oLthese results is obtained, 
by considering the more recent results from WUPE-3, LANLE3€3, UKQCD E3 
and APEEZI which all use either the improved SW action or the Wilson action 
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Figure 7: Lattice results for D — > K and D — * if* semileptonic decay form factors at zprp 
momentum transfer, and comparison with experimental results form the survey in Ref. Ej. 
The shaded bands show our summary values, indicating which results they are based on, as 

discussed in the text. 



with a KLM normalization. The summary values are presented in Table ^, 
along with recent experimental world averagesE3 in the D — > K, K* case. The 
values quoted reflect the fact that / + and A\ are the best measured while 
the D — > ir,p form factors are smaller with slightly larger errors. The shaded 
bands in Figure ^ show our estimates in the D — ► K, K* case and indicate 
which results they are based on. 

One sees that the lattice and experimental results agree rather well. The 
lattice values for A\ and A2 are both high compared to experiment: however, 
these depend on the correct normalization of the lattice axial vector current 
which is less well known than the vector current normalization needed for f + 
and V. In particular, the non-degeneracy of the c- and s-quarks means that 
there is no natural normalization condition to use for the weak current. This 
contrasts with the situation for heavy-to-heavy semileptonic decays, described 



in Section 4.2, where one can make use of the conservation of the vector current 



of degenerate quarks. 



4-2 Semileptonic B — > D and B — > D* Decays 

Semileptonic B — ► D* and, more recently, B — > D decays are used to determine 
the V c b element of the CKM matrix. Heavy quark symmetry is rather powerful 
in controlling the theoretical description of these heavy-to-heavy quark tran- 
sitions, as described by Neubert injlhe chapter on "B Decays and the Heavy 
Quark Expansion" in this volume O. In describing lattice results for j-these 
decays, we will quote theoretical results and direct the reader to Ref. O for 
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Table 6: Collected lattice results for D — > K , K* semileptonic decay form factors at q 



Yr /+(0) 



/°(0) 



V(0) 



V/Ai 



wue. 

LANLE3 
UKQCD 
APE 
ELC 
BKS 
LMMS 



97 

95-96 
95 
95 
94 

91-92 
89-92 



0.78(5) 

0.71(4) 

0.67(1) 

0.78(8) 

0.60(15)(7) 

0.90(8)(21) 

0.63(8) 



0.78(5) 
0.73(3) 
0.65(7) 



0.70(8)(24) 



1.27(16) 
1.28(7) 

i.oi(f§) 

1.08(22) 
0.86(24) 
1.43(45)(|i) 
0.86(10) 



1.78(7) 
1-4(1) 
1.6(3) 
1.3(2) 

1-99(22)® 
1.6(2) 



Yr Ai(0) A 2 (0) A 2 /A 1 A o (0) 



WIEE 
LANLB 


34 


97 




0.67(4) 


0.67(13) 




56 


95- 


96 


0.72(3) 


0.49(9) 


0.68(11) 0.84(3) 


UKQCD 


33 


95 




0-70^1) 


0.66(18) 


0.9(2) 0.75G5) 


APE 


37 


95 




0.67(11) 


0.49(34) 


0.7(4) 


ELC 


39 


94 




0.64(16) 


0.40(28)(4) 


0.6(3) 


BKS 


70 


91- 


92 


0.83(14)(28) 


0.59(14)(i) 


0.70(16)(?§) 0.71(16)(25) 


LMMS 


71 


89- 


92 


0.53(3) 


0.19(21) 


0.4(4) 



Table 7: Collected lattice results for D — > n,p semileptonic decay form factors at q 2 = 0. 



Yr /+(0) /°(0) V(Q) V/Ax 



WUP 


54 


97 


0.73(6) 


0.73(6) 


1.18(17) 




LANL 


35 


95 


0.56(8) 


0.62(5) 


1.18(15) 


1.77(16) 


UKQCD 


S3 


95 


0-61(1?) 


0.53(1?) 


0.95(ft) 




BKS 


70 


91-92 


0.84(12)(35) 


0.62(6)(34) 


1.07(49)(35) 


2.01(40)(32) 


LMMS 


71 


89-92 


0.58(9) 




0.78(12) 





Yr A x (0) A 2 (0) A 2 /A x A o (0) 



WUP 


34 


97 


0.63(5) 


0.64(14) 




LANL 


3.3 


95 


0.67(7) 


0.44(24) 0.67(31) 




UKQCD 


33 


95 


0.63(|) 


0.51(ig) 


0.70d|) 


BKS 


70 


91-92 


0.65(15)$) 


0.59(31)(I) 0.89(37)(?i) 


0.64(17)(21) 


LMMS 


71 


89-92 


0.45(4) 


0.02(26) 
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Table 8: Summary of lattice and experimental results for D — * K, K* and D — > it, p semilep- 
tonic decay form factors at q 2 = 0. The experimental numbers are taken from the survey in 

Ref.El. 





K,K* 


D — > 7T, p 




lattice 


expt 


lattice 


/ + (0) 


0.73(7) 


0.76(3) 


0.65(10) 


V(0) 


1.2(2) 


1.07(9) 


1.1(2) 


AM 


0.70(7) 


0.58(3) 


0.65(7) 


A 2 (0) 


0.6(1) 


0.41(5) 


0.55(10) 



details and references. 

In the heavy quark limit all six form factors in Eqs. (132-35) are related 
and there is just one universal form factor £(w), known as the Isgur-Wise (IW) 
function which contains all the non-perturbative QCD effects. Specifically: 



f+{q 2 ) = V(q 2 ) = A (q 2 ) = A 2 (q 2 ) 



(niB + mo) 2 



M{q 2 ) 



m B +m,D 



2^/mBrriD 



(37) 



where u> — vb-vd is the velocity transfer variable. Here the label D represents 
the D- or Z?*-meson as appropriate (pseudoscalar and vector mesons are de- 
generate in this leading approximation). Vector current conservation implies 
that the IW-function is normalized at zero recoil, i.e. that £(1) = 1. This 
property is particularly important in the extraction of V c b- 

The relations in Eq. ( |37|) are valid up to perturbative and power correc- 
tions. Allowing for corrections to the heavy quark limit, one writes the decay 
distribution for B — > D* as 



dr 

dui 



Q2 

-^(m B -ffl D .)V I) .\/w 2 -l(w+l) : 

4o7T J 



1 



4u> m 2 B — lujmBVtiB)- 4 
uj+1 (m B — m D ,y 



\V cb \ 2 T 2 (uj) 



(38) 



where !F(ui) is the "physical form factor" given by the IW-function combined 
with perturbative and power corrections. To extract \V c b\ one extrapolates 
measurements of the product ^(wJIV^tl to the zero-recoil point, w = 1, and 
uses a theoretical evaluation of the normalization :F(1)E3. 

A theoretical understanding of the shape of the physical form factor would 
be useful to guide the extrapolation of the experimental data, which currently 



26 



show rather a wide variation for the slope and interceptOQ, and also as a test 
of our understanding of the QCD effects. We expand J 7 as a power series in 
uj-V. 

T{uj) = [l - p 2 (u> - 1) + c (w - l) 2 + • • •] . (39) 

The slope parameter p 2 differs from the slope parameter p 2 of the IW function 
itself by heavy quark symmetry violating corrections!! 2 !, 

p 2 = p 2 + (0.16 ± 0.02) + power corrections. (40) 

To discuss lattice results for the shape of the IW function, and to search 
for corrections to the relations obtained using heavy quark symmetry, it is 
convenient to work with the following set of form factors (expressed in terms 
of four- velocities) which in the heavy quark limit either vanish or are equal to 
the IW function: 

(v+v'Th+(oj) + {v-v'Yh-{u) 

(oj+l)e*»h Al (<j) - e*-v(v»h A2 (u)+v'»h A3 (u)) 
where 

hi(u) = (an + AM + 7iM)£M (41) 

with 

a + = ay = a Al = a Aa = 1, Q<- = ua 2 = 0. (42) 

The (3i and 7, denote perturbative and power corrections (in l/m(, lC ) respec- 
tively. The statement of Luke's theoremO is 

7+ , Al (l)=0(A 2 QCD /< 6 ). (43) 

The principal difficulty for lattice calculations is to separate the physical 
heavy quark mass dependence due to power corrections from the unphysical 
one due to mass-dependent discretization errors. One must also address the 
question of lattice-to-continuum matching. We will illustrate our_discussion 
using the analysis procedure applied by the UKQCD collaborationLJ (who use 
the SW action for all quarks so that the leading mass-dependent discretization 
errors are formally reduced to 0(a s mQa) and 0(mgfl 2 )) . Consider the case of 
h + . For this form factor we have the protection of Luke's theorem at zero recoil 



(D(v')\cy>*b\B(v)) 







(D* 


{v')\c^b\B{v)) 






(D*(v' 


■)\cr 75 b\B(v)) 



Jm B m D 
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Table 9: Values of the Slope of the IW— function of a heavy meson, obtained using lattice 
QCD. The table indicates which functional form from Eq. Q47I) has been fitted to and which 



form factor has been used in the extraction. The systematic error in 
incorporates the variation from fitting to all functional forms in Eq. ( 



UKQCD results 
BSS note that 



is 12% smaller than p^, but do not quote a separate result. 



Yr 



LANL 
UKQCD 
UKQCD 
BSS 
BSS 



fit 



using 



96 0.97(6) 

95 0.9(|)(|) 

94 0.9(f) (?) 
93 
93 



NR h+ 

1.2(|)(?) NR h+ 

1.2(i) (I) NR h Al 

1.24(26)(36) lin h+ 

1.41(19)(41) NR h+ 



and, for degenerate (Q = Q') transitions, 
QlfiQ provides the further constraints: 



conservation of the vector current 



P+(l;mQ,m,Q)=Q, j + (l;mQ,m Q ) = 0. 



(44) 



The correct normalization of the vector current can be assured by requiring 
that the electric charge of the meson be 1. We therefore define the continuum 
form factor by, 



h + {ui\m Q ,mQ>) = (l + /?+(!; m Q, m Q')) 



h\(uj;m Q ,m Q ,) 
h\(l;m Q ,m Q >) 



(45) 



where h}"{ui; tuq, itiqi) is the (un-normalized) form factor calculated directly 
in the lattice simulation. This definition partially removes discretization er- 
rors and also removes w-independent power corrections while maintaining the 
known normalization conditions. If the remaining power corrections are small, 
then 



h+(u}) 



(46) 



1 + /3+M 

is effectively the IW function, This is convenient for extracting £(w), but 



the definition of Eq. ( |45[) precludes a determination of the zero-recoil power 
corrections. These corrections should be small, being suppressed by two pow- 
ers of the heavy quark mass. However, applying an analogous procedure to 
the Iia^lo) form factor relevant for B — > D* decays will not allow the 
corrections to .F(l), one of the dominant theoretical uncertainties, to be de- 
termined. | — | 

UKQCD confirmEa that their results for h+/(l + (3+) are indeed indepen- 
dent of the heavy quark masses and hence demonstrate, within the available 
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Figure 8: Fit of the UKQCD lattice results for |V.i.|J 7 ftiAff^ to the experimental data from 

the CLEO collaboration EE 



precision, that there is an IW function. -Moreover, a similar analysis for the 
h^! form factor reveals the same function^, so the IW function appears indeed 
to be universal. Performing extrapolations to the light (u-, d-) and strange (s-) 
quark masses and fitting to 

reNRH^^ T exp(-2(p 2 -l)^ T i) NR 
£(w) = < 1 - p 2 (uj - 1) linear (47) 

[ 1 — p 2 (uj — 1) + |(w - l) 2 quadratic 

gives lattice determinations of the slope of the IW function, as listed in Table ^[ 
Since 'the' IW function is different for different light degrees of freedom, the 
results in the table are labelled with subscripts u, d or s as appropriate- 
In Figure || we show the comparison of the UKQCEUlattice results l3 with 
B — > D* data (in 1995) from the CLEO collaborationlll. A chi-squared fit is 
made to the experimental data for |V^,|[1 + (3a ± (1)]-K"( w )£u,d( w ) = |^cfc|-?"(w). 
K(lu) incorporates radiative corrections away from zero-recoil (/3 J 4 1 (w)), non- 
perturbative power corrections (7^(0;)) and the contributions of the other 
form factors to the rate. The current lattice calculations cannot distinguish 
the u> dependence of K(oj), so it is taken to be a constant and hence p\ d and 
p 2 are not distinguished. The slope of the IW function is constrained to the 
lattice result in the fit so that the only free parameter is |Vci>|-^"(l). The result 
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of the fit is0 



\V cb \F{l) = 0.0371 \+l±\. 



(48) 



We should also mention B -^.EL semileptonic decays, which are beginning 
to be measured experimentally LJ'O with good precision, despite the helicity 
suppression in dT(B — > Dlvi)/du>. The differential decay rate depends on both 
h + {oj) and h_(ui). However, h_ is rather poorly determined to date in lattice 
calculations, so that it is difficult to evaluate the O(l/mo) corrections. 

Direct lattice calculations of the 1W function are-heing undertaken using 
discretizations of the heavy quark effective theory t3~E3. These studies are 
very interesting, but the results are not yet useful for phenomenology. An 
interesting theoretical feature of this approach is the formulation of the HQET 
at non-zero velocity in Euclidean space a~E3. 

Finally we note that a first lattice study of the semileptonic decays — * 
K c lv and — > H c lv has recently been performed ED, giving predictions for 
the decay distributions and the baryonic fsgur-Wise function. We refer the 
interested reader to Ref. EZl for details. 



4-3 Semileptonic B — > p and B — > 7r Decays and the Rare Decay B — ► K*j 

In this subsection we consider the heavy-to-light semileptonic decays B — > p 
and B — > 7r whieh are now being used experimentally to determine the V u b 
matrix element Several gxmr^Jiava. evaluated form factors fox these 

30~Eil (see the recent review inM) 



decays using lattice simulations liJ'Ej (see the recent review in U). We 

will also consider the rare radiative decay B — > K*j which is related by heavy 
quark and light flavour symmetries to the B — > p semilep. tqni p decay, and 
which was observed experimentally for the first time in I993t3t 



Form factors for semileptonic B — ► n decays were defined above in Eq. (33) 
and for B — > p decays in Eq. ([35]). For completeness, we define here form factors 
for the matrix element of the magnetic moment operator responsible for the 
short distance contribution to the B — > AT*7 decay: 

3 

(K*(k,e)\sa flv q v b R \B( P ))=J2ciT i (q 2 ), (49) 

i=l 

where q = p—k, e is the polarization vector of the K* and 

Cl = 2e^ A , £ *Vfc p , (50) 
Cl = £>B-»n]tO-e , ?(P + %. ( 51 ) 
Cl = e*- q L- q " ( P + k)X (52) 

\ ' B " l K* / 
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Table 10: Leading M dependence of form factors for heavy-to-light decays in the helicity ba- 
sis. The dependence follows from heavy quark symmetry applied at fixed velocity transfer u. 



Note that only three of the four A% form factors for B — > pli 


/ are independent. 


form i-channel leading M 


form 


i-channel leading M 


factor exchange dependence 


factor 


exchanj 


ie dependence 


B -> plv 




B -. 


■ -kIv 


V I- M 1 / 2 


/+ 


\- 


M 1 ' 2 


A x 1+ M- 1 ' 2 


/° 


0+ 


M- 1 ' 2 


A 2 1+ M 1 / 2 




B -> 




A 3 1+ M 3 / 2 


Ti 


1- 


AfVa 


A Q 0- M 1 ' 2 


T 2 


1 + 


M-Va 



T3 does not contribute to the physical B — > iiT*7 amplitude for which g 2 = 0, 
and Ti(0) and T 2 (0) are related by, 

T 1 (g 2 =0) = ^^ 2 ( (Z 2 =0). (53) 

Hence, for the process B — » -ftr*7, we need to determine Ti and/or T 2 at the 
on-shell point g 2 =0. 

Heavy quark symmetry is less predictive for heavy-to-light decays than 
for heavy-to-heavy ones. In particular, there is no normalization condition at 
zero recoil corresponding to the condition £(1) = 1, which is so useful in the 
extraction of V c b- The lack of such a condition puts a premium on the results 
from nonperturbative calculational techniques, such as lattice QCD. Heavy 
quark symmetry does, however, give useful scaling laws for the behaviour of 
the form factors with the mass of the heavy quark at fixed u>. -Moreover, the 
heavy quark spin symmetry relates the B — * V matrix elementsE-IES (where V 
is a light vector particle) of the weak current and magnetic moment operators, 
thereby relating the amplitudes for the two processes B° — > p + l~Di and B — > 
if *7, up to SU(3) flavour symmetry breaking effects. 

For fixed u> the scaling laws for the form factors given by heavy quark 
symmetry are as follows: 

Z(? 2 H) I fixcd e = M"/ 7/ (1 + % + ^ + • • •) (54) 

where / labels the form factor, M is the mass of the heavy-light meson and 
8 is a calculable leading logarithmic correction. The leading M dependences, 
M v f , are listed in Table Lattice calculations with propagating quarks use 
a range of quark masses around the charm mass and generally employ these 
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Figure 9: Ratios V/2T\ and Ai/2«T2 for five values of uj at three different heavy-light 
pseudoscalar masses, around the D mass, around the B mass and in the infinite mass limit. 
The horizontal dashed line denotes the heavy quark syrrmetry prediction in the infinite mass 

limit. Data from Ref. | 



scaling relations to extrapolate to the B mass: this is the case for results from 
ELC, APE and UKQCD. In the limit M — ► oo we also have the relations 

A x (q\u>)) = 2 l T 2 (q 2 (u J )), V(q 2 {oj)) = 2T 1 (g 2 H), (55) 

at fixed to. The UKQCD collaboration have checked the validity of the relations 
in Eq. (^5|) E3, finding that they are well satisfied in the infinite mass limit 
as shown in Figure ^. However, the ratio V/2T\ already shows significant 
deviations from the limiting value, 1, at the B mass. 

There are also kinematic constraints on the form factors at q 2 = 0: 

/+(o) = /°(o), r 1 (o) = »r a (o), a (o) = a 3 (o), (56) 

which will be useful below. 

From lattice simulations we can obtain the form factors only for part of 
the physical phase space. In order to control discretization errors we require 
that the three-momenta of the B, it and p mesons be small in lattice units. 
This implies that we determine the form factors at large values of momentum 
transfer q 2 — (ps — Pir,p) 2 ■ Experiments can already reconstruct exclusive 
semileptonic b — > u decays (see, for example, the review inl23) and as techniques 
improve and new facilities begin operation, we can expect to be able to compare 
the lattice form factor calculations directly with experimental data at large q 2 . 
A proposal in this direction was made by UKQCDEJ for B° — > p + l~Di decays. 
To get some idea of the precision that might be reached, they parametrize the 
differential decay rate distribution near g^ax by: 

^2 = 10 T09^3Tp 1 X2 (l ) a { 1 + b (l - W)) . ( 57 ) 
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q 2 (GeV 2 ) 

Figure 10: Differential decay rate as a function of q 2 for the semileptonic deca y B ° — > 
p+i — P;, taken fromEj. Points are measured lattice data, solid curve is fit from Eq. ( ply with 
parameters given in Eq. (p^). The dashed curves show the variation from the statistical 
errors in the fit parameters. The vertical dotted line marks the charm endpoint. 



where a and b are parameters, and the phase-space factor A is given by X(q 2 ) = 
{m 2 B + m 2 — q 2 ) 2 — Arn^rn 2 . The constant a plays the role of the IW function 
evaluated at u> = 1 for heavy-to-heavy transitions, but in this case there is no 
symmetry to determine its value at leading order in the heavy quark effective 
theory. UKQCD obtain E 

a = 4.6±8:t±0.6GeV, 
b = (-8+f) x l(r 2 GeV 2 . 

The fits are less sensitive to b, so it is less well-determined. The result for 
a incorporates a systematic error dominated by the uncertainty ascribed to 
discretization errors and would lead to an extraction of \V u b\ with less than 
10% statistical error and about 12% systematic error from the theoretical input. 
The prediction for the dT/dq 2 distribution based on these numbers is presented 
in Figure [HJ With sufficient experimental data an accurate lattice result at a 
single value of q 2 would be sufficient to fix \V u b\- 

In principle, a similar analysis could be applied to the decay B° — * Tt + l~T>i. 
However, UKQCD find that the difficulty of performing the chiral extrapolation 
to a realistically light pion from the unphysical pions used in the simulations 
makes the results less certain. The B — > n decay also has a smaller fraction of 
events at high q 2 , so it will be more difficult experimentally to extract sufficient 
data in this region for a detailed comparison. 

We would also like to know the full q 2 dependence of the form factors, 
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Figure 11: Bounds on / + ar>d-i/° for B° — » 7r+£~P; from dispersive constraints S. The 
data points are from UKQCDEJ, with added systematic errors. The pairs of fine curves are, 
outermost to innermost, 95%, 70% and 30% b«unds. The upper and lower shaded curves 
are light-cone E3 and three-point E3 sum rule results respectively. 



which involves a large extrapolation in q 2 from the high values where lattice 
calculations produce results, down to q 2 = 0. In particular the radiative decay 
B — > K*j occurs at q 2 = 0, so that existing lattice simulations cannot make a 
direct calculation of the necessary form factors. 

An interesting approach to this extrapolation problem has been applied 
by LellouchE3 for B° — > ir + l~9i. Using dispersion relations constrained by 
UKQCD lattice results at large values of q 2 and kinematical constraints at q 2 = 
0, one can tighten the bounds on form factors at all values of q 2 . This technique 
relies on perturbative QCD in evaluating one side of the dispersion relations, 
together with general properties of field theory, such as unitarity, analyticity 
and crossing. It provides model-independent results which are illustrated in 



Figure 11. The results (at 50% CL — see Ref. for details) are 



/+(0) = 0.10-0.57, (59) 
T(B° -> 7r+r h) = 4.4-13 K 2 &|PS _1 . (60) 

Unfortunately these bounds are not very restrictive when constrained by exist- 
ing lattice data. In principle, this method can be applied to B — > p decays also, 
but is more complicated there, and the calculations have yet to be performed. 
Recently, Becirevicl£3 has applied the method for B — > K*~/, using APE lattice 
results as constraints. However, he has not applied the kinematic constraint 
from Eq. (Bq) and the resulting bounds are not informative: they become so, 
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Table 11: Lattice results for B° —> 
The decay rates are values for T(B° 



7r+i u t using various ansatzp-4br the for* 
r+J-piVIVuspps- 1 . ELCE2 and APEE 



I factor /+. 
I results are 



from their method 'b', which uses the heavy quark scaling laws of Eq. (m) to extrapolate 

from D- to B-mesons at fixed lo. 



Yr (3 



norm Rate / + (0) 



UKQCD 
WUP 
APE 
ELC 



97 6.2 1 rel 

97 6.3 nr 

95 6.0 1 rel 

94 6.4 rel 



8.5(ft) 0.27(11) 
0.43(19) 
8±4 0.35(8) 
9± 6 0.30(14)(5) 



however, once he uses a light-cone sum rule evaluation of Ti(0) = 2X2(0) as an 
additional constraint. These dispersive methods can be used with other ap- 
proaches in addition to lattice results and sum rules, such as quark models, or 
even in direct comparisons with experimental data, to check for compatibility 
with QCD and to extend the range of results. 

For now we must rely on model input to guide q 2 extrapolations. We can 
ensure that any assumed g 2 -dependence of the form factors is consistent with 
the requirements imposed by heavy quark symmetry, as shown in Eq. (p4|) , to- 
gether with the kinematical relations of Eq. (pfj). Even with these constraints, 
however, current lattice data do not by themselves distinguish a preferred q 2 - 
dependence. Fortunately, more guidance is available from light-cone sum rule 



analyscsl \\ ^ which lead to scaling laws for the form factors at fixed (low) q 



2 



rather than at fixed u as in Eq. ([>4|). In particular all form factors scale like 
M- 3 / 2 at q 2 = 0: 



/(0 ,e = M-^ 7 ,(i+^ + ^ + 



(61) 



It is therefore important to use ansatze for the form factors compatible with 
as many of the known constraints as possible. 

Lattice results for B° — > tt + 1~Di, B° — ► gd.l~vi and B\-^r> K*l are reported 
in Tables [ll], |l| and |l| respectively. ELC@ and APE ED fit lattice data for 
the semileptonic decays at a single value of q 2 to a simple pole form with the 
appropriate pole mass also determined by their data. For the f° and A\ form 
factors, this is consistent with heavy quark symmetry requirements, kinematic 
relations and light-cone scaling relations at q 2 = 0, Eq. ([mLbut for the other 
form factors it is not simultaneously consistent. The WUPEj results are found 
by scaling form factors at q 2 — from results with quark masses around the 
charm mass to the 6-quark mass. However, the scaling laws used do not follow 
the light-cone scaling relations of Eq. mW) . 
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Table 12: B° — ► p+i ^; results frpwi lattice simulations. The decay rates are values for 
T(B° -> p+i— S>j) / 1 1 2 ps~ 1 . ELCH andAPEED results are from their method 'b', which 
uses the heavy quark scaling laws of Eq. (m) to extrapolate from D- to B-mesons at fixed u>. 







Yr 






norm 


Rate 


V(0) 


Mo) 


A 2 (0) 


UKQCD 
WUP 
APE 
ELC 


n 

34 
37 
30 


97 
97 
95 
94 


6.2 
6.3 
6.0 
6.4 


1 



1 




rel 
nr 
rel 
rel 


16.5(1) 

12 ±6 
14 ±12 


0.35(1) 
0.65(15) 
0.53(31) 
0.37(11) 


0-27(1) 
0.28(3) 
0.24(12) 
0.22(5) 


0.26(|) 
0.46(23) 
0.27(80) 
0.49(21)(5) 



The latest UKQCD study Ell uses models consistent with all constraints, 
including the light-cone sum rule scaling relations^]. For B° — >■ Tr + l~i?i UKQCD 
use a dipole(pole) model for / + (/°), while for B° — ► p + l~vi and B p- K*j 
(collectively denoted by B — > V) they use a model inspired by Stecht£3 which 
expresses all the form factors for the scmilcptonic and radiative decays of a 
heavy pseudoscalar meson to a light vector meson in terms of a single function, 
taken to be a simple pole form for A\. In lattice calculations one has the 
freedom to adjust hadron masses by tuning the quark masses used in the 
simulation. UKQCD use this freedom to perform a combined fit within their 
model for all the B — > V form factors simultaneously, first with V — p and then 
with V = K* . From the first fit they obtain the form factors for B° — > p + l~vi 
and from the second they obtain results for B — > K*^. The combined fit in 
the K* case is illustrated in Figure |lj. The figure demonstrates the large 
extrapolation needed to reach q 2 = 0. 

Our preferred results-for B° — * n + l~vi and B Q — > p + l~Di come from the 
UKQCD constrained fitsEJ. Their values for the form factors extrapolated to 
q 2 = agree well with light-cone sum rule calculations, which work best at low 
q 2 . The fitted form factors also agree with experimental results for the rates 
and ratio-of-rates of these semileptonic decays. However, we emphasise that 
the extrapolated form factors are no longer model independent. 

There are also preliminary results for heavy-to-light form_£aL±ors from 
FNAL, JLQCD and a Hiroshima-KEK group (see the reviews inE2nij) and the 
different lattice calculations are in agreement for the form factors at large q 2 
where they are measured. 



Table |13| lists the values of T(0) = Ti(0) = iT 2 {0) for B -> K*j, together 
with the directly measured T^max)- All groups find that T 2 has much less 



*In the extrapolation of the raw lattice results to the B-meson mass, this analysis also 
applies the condition in Eq. that the pairs (V, 2Ti) and (A\,2iT2) should agree at 

fixed u) in the infinite heavy mass limit. 
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q 2 dependence than Ti, although the lattice data again do not themselves 
distinguish a preferred q 2 dependence. In order to make a distinction, one can 
apply the light-cone sum rule scaling relation at q 2 — 0, see Eq. (^TJ), which 
states that T(0) has a leading M -3 / 2 behaviour. In the table we list results 
from form factor fits which satisfy this scaling law. Our preference is to quote 
the UKQCD result (with statistical error only) from the combined fit to 
B —f V decays described above: 



T(0) 



0.16(f). 



(62) 



Using this value to eval 



to 0(l/mi) corrections EHrJ) 



e the ratio (given at leading order in QCD and up 



Rk' 



r(B^K* 7 ) _ / m L ; 



ni b 



\T(0)\< 



(63) 



results in 

R K . = 16(|)%, (64) 

which is consistent with the experimental result 18(7)% from CLEO@. Dis- 
crepancies between Rk* calculated using T(0) and the experimental ratio 
T(B — > K*j)/T(b —> sj) could reveal the existence of long-distance effects in 
the exclusive decay. It ha^ijeftuproposed that these effects may be significant 



for the process B — > K*j 



loe 



but within the precision of the experimental 



and lattice results, there is no evidence for them. 
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Table 13: Lattice results for B -> K*^. Values for T(0) = 71(0) = iT 2 (Olaxe quoted only 
from models which satisfy the light-cone sum rule scaling relation, Eq. (pll), at q 2 



0. 



Yr P 



norm 



IjkqcdE 



r(o) 



LANL 
APE 
BHS 



107 
108 
109 



97 6.2 1 rel 

96 6.0 nr 

96 6.0 1 rel 

94 6.0 nr 



0.16(f) 0.25(2) 

0.09(1) 

0.09(1)(1) 

0.101(10)(28) 0.325(33)(65) 



5 The Parameters of the HQET 

The Heavy Quark Effective Theory (HQET) is proving to be a particularly 
useful tool for phcnomenological studies in charrn.-a.nd .beauty physics, for re- 



views and references to the original literature seeOti3'tiil. In this approach, 
physical quantities are calculated as series in inverse powers of the mass(es) 
of the heavy quark(s). The non-perturbative strong interaction effects can be 
parametrized in terms of matrix elements of local operators, which appear as 
factors in the expansion coefficients. Lattice simulations of the HQET provide 
the opportunity of computing these matrix elements numerically and in this 
section we briefly describe some of these calculations. We start with a gen- 
eral discussion (for_a-more detailed presentation and references to the original 



literature see Ref. 1113) and then we will illustrate the main points by some 
important examples. 

Consider some physical quantity V which, using the operator product ex- 
pansion can be written as a series in inverse powers of the mass of the heavy 
quark, m Q , |] 

V = ^(my^X/IOi W) + C2(m ^ 2) (/l<WN) + 0(l/m n Q +p ) , (65) 

Q 

where n > 1 and the coefficient functions C, are independent of the states |z) 
and |/). The operators 0\.i are composite operators of static heavy quark, 
light quark and gluon fields, renormalized at a scale /i. For clarity of nota- 
tion, we suppress the dependence of the coefficient functions on the coupling 
constant, a s (rriQ). We assume here that there is only one operator in each of 
the first two orders of the expansion. If this is not the case, then there is an 
additional mixing of operators, which requires only a minor modification of the 
discussion below. We will therefore not consider this possibility further. The 



s For simplicity we assume that V depends only on the mass of one heavy quark. 
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final term of 0(l/mg +p ) in Eq. ( |65| ) represents the contributions of operators 
of even higher dimension which will not be discussed here. 

Throughout this discussion we assume that we wish to evaluate the con- 
tributions of 0( (Aqcd/toq)") to V, or at least to improve the precision given 
by the lowest order contribution. 

In lattice simulations we work directly with bare operators (corresponding 



to the lattice regularization), so that the renormalization scale [i in Eq. (65) 
should be replaced by the cut-off a -1 . With a hard ultra-violet cut-off, such 
as the lattice spacing, the higher dimensional operators 02,3,... can mix with 
lower dimensional ones, and with 0\ in particular. By dimensional arguments 
it can readily be seen that the mixing coefficients will diverge as inverse powers 
of the lattice spacing. This implies, for example, that 

(f\0 2 (a)\i)=0(a- n ) , (66) 

and so this matrix element, computed in lattice simulations, is clearly not a 
physical quantity. These power divergences are subtracted in the perturbative 
matching procedure in which the Wilson coefficient functions are computed; 
specifically the coefficient function C\ contains terms of 0((ma)~ n ), which 
cancel the power divergences present in (/| O2 The perturbative series for 
C\ is evaluated only at some low order (typically one loop). Since (1/mQd)™ is 
much larger than (Aqce> / m o) n , which is the order of the terms we are trying to 
evaluate, the truncation of the perturbative series for C\ leads to a significant 
loss of precision This is the principal problem, in attempts to quantify 
power corrections to hard scattering and decay amplitudes in general. We 
mention in passing that although the discussion here is in the context of lattice 
calculations, analogous problems occur also_in the continuum regularization 
schemes, such as dimensional regularizationt 
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In some cases the operator O2 is protected from mixing under renormal- 
ization with 0\ , because of the presence of some symmetry. An important ex- 
ample of this in heavy quark physics is the chromomagnetic operator haijG lJ h 
(where h is the field of a static quark and G^ v is the gluon field strength ten- 
sor), which cannot mix with the lower dimensional operator hh because it has 
a different spin structure. In such cases the corresponding problem of power 
divergences (and renormalon singularities) does not arise. For the remainder 
of the discussion we assume that O2 and 0\ have the same quantum numbers, 
so that they can mix under renormalization. Another important exception to 
the general discussion is the difference of matrix elements of the kinetic energy 
operator taken between different hadronic states. In this case the higher di- 
mensional operator, KD 2 h, can mix with the lower dimensional one, hh, but 

-0. 



''The connection of the power divergences and renormalons is explained in Refs. 
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as the latter is a conserved current it has the same matrix element between all 
single hadron states. Thus the corresponding power divergences cancel in the 
difference of matrix elements. In all of these exceptions the matrix elements of 
the higher dimensional operators (or linear combinations of matrix elements) 
are also the leading contribution to some physical quantity (for example, the 
chromomagnetic operator gives the leading contribution to the B*-B mass 
splitting). 

In the following subsections we consider 3 parameters which occur fre- 
quently in the evaluation of physical quantities using the HQET: 

• A, the binding energy of a heavy quark in a heavy hadron; 

• Ai , the kinetic energy of the heavy quark; 

• A2, the matrix element of the chromomagnetic operator. 

The binding energy A is defined as the difference of the masses of the heavy 
hadron and the heavy quark. In practice, the heavy quark mass must be defined 
at short distances (otherwise experimentally measurable quantities cannot be 
expressed in terms of the mass using perturbation theory), and below we will 
use the MS mass. The usefulness of introducing A is then unclear, and we 
present the discussion simply in terms of the mass itself. In all the three 
examples to be discussed below, results are presented from lattice calculations 
using static heavy quarks. 



5. 1 The Evaluation of the Mass of a Heavy Quark 

The first example which we consider is the determination of the mass of a 
heavy quark (e.g. the 6-quark), up to, and including the terms of O(Aqcd), 
but neglecting terms of 0(Aqq D / m^) j|. 

In section ^, the computation of the decay constant of a meson containing 
a static heavy quark was discussed. This parameter is obtained by evaluating 
correlation functions of the form: 

C(t) = ^<0|A 4 (x,t)A4(0,0)|0) (67) 

X 

in lattice simulations of the HQET, with Lagrangian density 

C = hD A h , (68) 



*A formulattmirofithis problem in terms of an explicit operator product expansion can be 
found in Refs.£I3EJ. 
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where h represents the field of the static quark. For sufficiently large values of 
the time t, 

C{t)~Z 2 e- £t + --- , (69) 

where the ellipsis represents contributions from excited states. The value of 
fs in the static limit is obtained from the prefactor Z. In addition, however, 
from the exponent £ it is possible to obtain the O(Aqcd) contribution to m^. 
Performing the matching of the heavy quark propagator in full QCD and in 
the HQET givesN 



£ = m b — ( m b° lc ~ & m ) t (70) 



where raj is the mass of the B-meson and m^ olc is the pole mass of the b- 
quark. Although we have chosen the bare Langrangian ( |68| ) to have no mass 
term, higher order perturbative corrections generate such a term and 

n=l 

represents the perturbation series generating the mass. Both £ and 8m di- 
verge linearly with the lattice spacing, so that they are not physical quantities. 

nnphysical, and in perturbation theory contains renormalon am- 



polc 

mr b is qJSQ 



biguities t L15 H 117 l These ambiguities cancel those also present in 8m in the 
difference on the right-hand side of Eq. ([70|) . Thus it is possible to determine 
the value of a physical (short-distance) definition of the quark mass, such as 
To;, = m^ s (m^ IS ) from the computed value of £. In practice, however, the 
subtraction of the linear divergence, which is performed in perturbation theory 
leads to larg£_numerical cancellations and hence to significant uncertainties. 
In reference tlB it was found that 



m b = 4.15 ±0.05 ±0.20 GeV. (72) 

The first error on the right-hand side of Eq. ( [72] ) is due to uncertainties in 
the lattice evaluation of £ and in the value of the lattice spacing. However, 
the larger error of 200 MeV or so is the estimate of the uncertainty due to 
the truncation of the perturbation series for 8m at one-loop order. Evaluation 
of higher_.order terms in this series, perhaps using the methods developed 
based on the Langevin stochastic formulation of lattice QCD, is 



m 



Ref. 



urgently needed to reduce the uncertainty in the computed value of m^. 
5.2 Kinetic Energy of a Heavy Quark 

The next example which we will consider is the evaluation of the kinetic energy 
of the heavy quark Ai, which appears in many applications of the HQETL3. 
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On the lattice we start with the evaluation of the matrix element of the bare 
operator 

A?" = - (73) 
1Mb 

= -(0.69 ± 0.03 ± 0.03) a- 2 (74) 

where the numerical result is taken from a simulation on a 24 3 x 60 lattice 
at P = 6.0 with the SW-action Taking the lattice spacing to be 2.0 ± 
0.2 GeV, we see that the magnitude of the result is about 2.8 GeV 2 , to be 
compared to the expected physical corrections of 0(Aq CD ) ~ 0.1 GeV 2 . Of 
course, the large result is due to the presence of power divergences, in this case 
they are quadratic, i.e. they are of 0(a~ 2 ). In one-loop perturbation theory, 
the power divergence is equal to — 5.19 a s a -2 , which, depending on the value 
taken for the coupling constant a s , is in the range (0.67 — 0.93)a -2 (the choice 
of a suitable definition of the coupling constant is a representation of our 
ignorance of the higher order perturbative corrections, the range given here 
comes from frequently used definitions). The uncertainty is greater than the 
terms we are trying to evaluate which are of 0(Aq CD ). Clearly, in order for the 
lattice result to be useful for phenomenological applications, the perturbative 
calculations must be performed to higher orders, which is the main conclusion 
of this subsection. 
I 

Refs. 



isalso ppssiblp to subtract the power divergences non-perturbatively. In 
a subtracted kinetic energy operator, 



KDlh=KD 2 h- Arhh , (75) 



was defined, with the subtraction contant c fixed by imposing that the matrix 
element of this operator vanishes between quark states at rest (in the Lan- 
dau Gauge). The corresponding value of Ai, which is now free of quadratic 
divergences, was found to be-[] 

Ai = a _a Z D |(a a Aj are -a a c) (76) 
= 0.09 ±0.14 GeV 2 . (77) 



Of course the large relative error in Eq. (|77|) is due to the large cancelation 
between the two terms in the parentheses in Eq. (f76|). In Eq. (|76|) , Z D 2 is the 
normalization constant required to obtain the continuum, MS, value of Ai from 
the subtracted lattice one. 



J Note that the central value in Eq. ( F^ ) has the opposite sign to that of many other 
estimates using various definitions of Ai. 
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The difficulties described here arise because of the mixing of the kinetic 
energy operator h(iT)) 2 h with hh/a 2 . Since hh is a conserved current in the 
HQET, its matrix elements are the same between all hadronic states. This 
means that the difference of the matrix elements of the kinetic energy opera- 
tor between any two different beauty hadrons is a physical quantity, and for 



example in Ref. E£2I it was found that 

Xi(B s ) - Xx(B d ) = -0.09 ± 0.04 GeV 2 . (78) 

This difference is the leading contribution to the following combination of 
masses 

X 1 (B S ) - \\{B d ) = m B m -ms d -(m Dt -m Dd ) + Q ( A|cd\ (7g) 







2 \ ran 


TUB J 



where, for example, is the spin averaged mass of the B- meson (with the 
corresponding light valence quark) 

m B = i(3m B » + m B ) , (80) 

and m,Q is the mass of the heavy quark Q (Q — b or c) . Where appropriate, we 
have included the subscript d or s to denote the presence of the corresponding 
light valence quark. The experimental value of the first term on the right hand 
side of Eq. (|7^) is —0.06 ±0.02 GeV 2 , in very good agreement with the result in 
Eq. (jr^) . It must however be remembered that the 0(AQ CD /m c ) corrections 
to this result may be significant. 

5.3 A2; The Matrix Element of the Chromomagnetic Operator 

The chromomagnetic operator h^aijG^ h, does not mix with the operator hh, 
and hence its matrix elements are free of power divergences. The parameter 
A2, defined as 

x ^^l^ h r« Gijh \ B) (81) 

gives the first term in the hyperfmc splitting in the i?-meson system: 

m|, - m| = 4A 2 . (82) 

The lattice results for A2 are all significantly smaller (by almost a factor of two) 
than the values deduced from the physical masses of the B* and B mesons. In 
two recent lattice computations the authors found: 



0.28 ±0.02 ±0.04 GeV 2 Ref. 54 
0.28 ±0.06 GeV 2 Ref. 



120 



(83) 



43 



to be compared to the experimental value of 0.485 ± 0.005 GeV 2 , see Ref.li^. 

One possible source for the discrepency between the lattice results and 
the experimental value is the unusually large one-loop contribution to the 
renormalization constant relating the lattice and continuum chromomagnetic 
operator E£3. This renormalization constant is about 1.85 at one-loop order, 
and so one may wonder whether the higher order terms might give a significant 
contribution. Other possible contributions to the discrepency might be the 
use of the quenched approximation, or that the relation between A2 and the 
hyperfine splitting may be significantly modified by higher order corrections in 
l/rub- It is important to clarify the source of this discrepency. 

Lattice calculations of the hyperfine splitting using propagating heavy 
quarks also give a result which is smaller than the experimental one. This 
is a different problem, however, which is related to the presence of a spurious 
chromomagnetic term of 0(a) present in the lattice action. This interpretation 
is confirmed by the fact that the computed value of the splitting increases as 
the action is "improved" in agreement with expectations I 123 !! 124 !. 



6 Exclusive Non-Leptonic Decays of Heavy Mesons 

Exclusive non-leptonic decays are, in principle, an important source of funda- 
mental information on the properties of weak decays of heavy quarks. Unfor- 
tunately our current theoretical understanding of the non-perturbative QCD 
effects in these processes is rather primitive, and we ai^ forced to make assump- 

ttice computations 



.1.2:, 



tions based on factorization and/or quark models 
of the corresponding matrix elements are also difficult They need to be 
performed in Euclidean space, where there is no distinction between in- and 
owi-states. The quantities which one obtains directly in lattice computations 
are the (real) averages such as 

(M 1 M 2 \H W \B) = ^( in ( Mi M 2 \H W \B}+ out < M x M 2 \ H w \B ) ) , (84) 

where Mi and M 2 are mesons. It is therefore not possible to obtain directly 
any information about the phase due to final state interact ions , and hence 
to determine the matrix elements reliably. Maiani and Testa E££l also showed 
that the quantities which are obtained from the large time behaviour of the 
corresponding correlation functions are the unphysical form factors in which 
the final state mesons are at rest, e.g. 

( Mi(p Ml = 0) M 2 { P m 2 = 0) \ Hw \B) . (85) 
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For K — > 7T7T decays chiral perturbation theory can then be used to obtain the 
physical form factors with reasonable accuracy. For D- and £?-meson decays 
this is not possible. 



The publication of the Maiani- Testa 



L26 



theorem stopped the exploratory 
work on the numerical evaluation of two-body non-leptonic decays. These 
early, and not very accurate-papers, studied the non-penguin contributions to 



127 



D-meson decay amplitudes 

The Maiani- Testa theorem implies that it is not possible to obtain the 
phase of the final state interactions without some assumptions about the am- 
plitudes. The importance of developing reliable quantitative techniques for 
the evaluation of non-perturbative QCD effects in non-leptonic decays cannot 
be overstated, and so attempts to introduce "reasonable" assumptions to en- 



able calculations to be performed-L 



compared with experimental data) are 
needed urgently. Ciuchini et al. tii3>E£if| have recently shown that by making a 
"smoothness" hypothesis about the decay amplitudes it is possible to extract 
information about the phase of two-body non-leptonic amplitudes. Studies to 
see whether their proposals are practicable and consistent are currently begin- 
ning. 



7 Conclusions 

The decays of heavy quarks, and the &-quark in particular, provide a power- 
ful laboratory to test the standard model of particle physics, to determine its 
parameters (in particular, the elements of the CKM-matrix) and to look for 
signatures of new physics. Non-perturbative QCD effects in these weak decays 
are the major difficulty in interpreting the present and future experimental 
data. Lattice computations provide the opportunity to quantify these effects, 
and much effort is being devoted to the evaluation of the decay amplitudes 
of heavy quarks. In this article we have reviewed the status of these calcula- 
tions, presenting an overview of the results for leptonic decay constants, the 
amplitudes for B — B mixing, the form factors for semileptonic and rare radia- 
tive decays and the parameters A, Ai and A2 of the HQET. These successful 
calculations are playing a central role in phenomcnological studies. 

As explained in Section ^, relatively little progress has been achieved so far 
in understanding exclusive non-leptonic decays on the lattice, or indeed any 
physical process in which there are two hadrons in the initial or final state. It 
is very important to try to overcome the theoretical difficulties which have so 
far prevented this very important class of processes from being amenable to 
study in lattice simulations. 

Lattice simulations are a "first principles" non-perturbative technique for 
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the evaluation of the strong interaction effects in weak decay amplitudes (and 
in those of heavy quarks in particular). However, the precision of the results 
is limited by the available computing resources, and the priority is to reduce 
and control the systematic uncertainties, especially those d ue t o the quenched 
approximation. This is discussed in some detail in Section 1.2 above. During 
the next few years, studies with dynamical fermions will continue and improve, 
enabling us to begin seriously quantifying the effects of the quenched approx- 
imation. The very significant progress achieved in recent years in reducing 
the errors due to the finitencss of the lattice spacing, makes this task eas- 
ier, allowing for meaningful studies on coarser lattices than would otherwise 
be the case. Once precise studies with dynamical quarks become possible, 
lattice QCD together with large scale numerical simulations will be a truly 
and fully quantitative technique for the evaluation of the effects of the strong 
interactions, with no model assumptions or parameters. Although this will 
still take a few years, we have tried to demonstrate in this review that lattice 
simulations, with their expected 0(10-20%) systematic uncertainties for most 
physical quantities, are already a most reliable and useful tool. 
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